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Summary 

A Controlled Ecological Life Support System (CELSS) 
utilizes a plant s natural ability to regenerate air and water 
while being grown as a food source in a closed life 
support system. Current plant research is directed toward 
obtaining quantitative empirical data on the regenerative 
ability of each species of plant and the system volume and 
power requirements. Two techniques were adapted to 
optimize crop species selection while at the same time 
minimizing the system volume and power requirements. 
Each allows the level of life support supplied by the plants 
to be selected, as well as other system parameters. The 
first technique uses decision analysis in the form of a 
spreadsheet. The second method, which is used as a 
comparison with and validation of the first, utilizes 
standard design optimization techniques. Simple models 
of plant processes are used in the development of these 
methods. 

Introduction 

To date, life support technology is based solely on 
physical/chemical processes, and this is likely to remain 
true for the initial phases of the Space Exploration 
Initiative. However, for long-duration missions, such as a 
trip to Mars or long-term habitats on the Moon or Mars, a 
Controlled Ecological Life Support System (CELSS) has 
the potential to provide human life support with signifi- 
cant cost and safety benefits over the currently envisioned 
physical/chemical systems. In particular, the amount of 
food resupplied from Earth may be significantly dimin- 
ished, higher plants can accomplish both air revitalization 
(through the release of oxygen and uptake of carbon 
dioxide) and water processing (through transpiration), and 
some waste disposal may be accomplished biologically. 
Figure 1 shows an example of an integrated biological and 
physical/chemical life support system. Other studies have 
resulted in variations on this conceptual design (refs. 1-3). 

Normally, a trade study is conducted to determine 
advantages and disadvantages of various design options. 
Trade-study techniques can be developed in parallel with 
research on basic performance parameters, so that when 
reliable data become available, the analysis tool is also 
ready to perform trades. These tools will become increas- 
ingly important as we begin to address the complexities 
involved in integrating biological components with 
physical/chemical life support system components. 

To date, most research in the use of plants for life support 
has concentrated on productivity levels and the effects of 
environmental parameters on productivity. Little work has 
been done in evaluating the air and water regeneration and 
waste management capabilities, which would be the next 


logical step toward developing an integrated life support 
system. The techniques outlined herein transform newly 
acquired plant performance data into parameters describ- 
ing a CELSS for use in trade studies, thus providing the 
link between generating data and developing an optimal 
CELSS design. 

We are adapting two techniques to optimize crop selection 
for minimum power and volume penalties. The first tech- 
nique involves the use of decision analysis which imple- 
ments a decision tree. The second method, intended to be 
both a comparison with and a validation of the first, 
involves the use of standard design optimization (linear 
programming) techniques. Previously, design optimiza- 
tion techniques have been applied to crop mix selection 
to select the minimum crop area which satisfies human 
nutritional requirements (ref. 4). While the results of that 
study do not account for power penalties, nor do they 
allow for air and water regeneration constraints, some 
comparisons can be made with our own data. These 
comparisons are reported in the Results section of this 
paper. 

A spreadsheet is used as an interface with the user and to 
generate plant and system parameters. The user specifies 
the level of life support to be supplied by the plants for 
each life support function— oxygen generation, carbon 
dioxide uptake, water regeneration, and nutrient 
(carbohydrate, lipid, and protein) production. Plant 
parameters are generated both directly from empirical 
data input by the user and from models of plant processes. 
However, the techniques are intended to be generic and 
applicable to other plant parameter generation schemes. 

Approach 

Derivation of Plant and System Parameters for Input 
to Analyses 

Parameters specifying system requirements are input into 
the spreadsheet section shown in table 1. Here the user 
can input which crops of those whose performance data 
have been entered into table 2 should be considered. Up to 
five crops can be chosen. If variety in the crop mix is 
desired, a minimum number of crops can be entered (a 
number greater than one will force multiple crops to be 
chosen). Currently this feature applies only to the decision 
analysis method. The level of life support to be fulfilled 
by plants and the daily life support requirements (ref. 5) 
are input into table 1, as well as the power and volume 
penalties (refs. 6 and 7) and end-to-end lighting 
efficiency. 



Table 2 lists parameters which were obtained from the 
literature and other sources of plant data (refs. 8-12). 1 
These parameters include each species * rate of tran- 
spiration and biomass growth as well as diet composition, 
edible fraction, and chamber height requirement. The 
lighting level and photoperiod under which these rates 
have been measured are also recorded and used for the 
power requirement calculations. There are a myriad of 
other factors that influence plant productivity, tran- 
spiration rates, edible fraction, and even diet composition, 
such as carbon dioxide level, oxygen level, nutrient 
solution composition, temperature, and humidity. These 
influences could be added in a more sophisticated effort to 
derive the parameters for the analyses input, but are not 
necessary for our purposes of technique development and 
demonstration. 

Table 3 lists the parameters describing plant species' 
performance required by the optimization methods. 
Generation of these parameters can be accomplished in 
many different ways, from using empirical data to 
employing modeling techniques. We have elected to use a 
combination of these two methods, largely because gas 
exchange data for plant species are limited. Parameters 
which are more readily available from empirical data are 
used directly and as the basis for some simplified relations 
used to generate the remaining parameters. 

Transpiration rate and biomass production rate are taken 
directly from the empirical data recorded in table 2. The 
fat, carbohydrate, and protein production rates are prod- 
ucts of the biomass generation rate, m b i 0 , the edible 
fraction, and the fraction of the total biomass generation, 
which is fat, carbohydrate, and protein, respectively. 


fat 


fat 

generation 


fraction 

carbohydrate 

generation 

+ ( edible 3 * 
” m bio fraction J 

carbohydrate 

fraction 

protein 

generation 


protein 

fraction 


Carbon dioxide and oxygen generation rates are products 
and reactants of photosynthesis and respiration, as is 
biomass production. This link between biomass 
production and gas exchange rates is described by the 
photosynthetic equation, assuming respiration is ignored. 


! Crops X and Y in table 2 have been added to enlarge the 
database. These two “species” are not real; the data from them 
merely illustrate the development of this tool. 


The chemical reaction of photosynthesis varies with the 
biomass type being formed, whether it is fat, carbo- 
hydrate, or protein. For simplicity, here it is assumed that 
carbohydrate is the substance formed. The photosynthetic 
equation describing carbohydrate formation is (ref. 13) 

6CO2 + 6H2O — » C5H12O6 + 6O2 

This equation gives the ratio of moles of carbohydrate 
produced to moles of carbon dioxide taken up and oxygen 
released. Converting these mole fractions to mass 
fractions (using the corresponding molecular weights), the 
oxygen generation rate, riiQ 0 , and the carbon dioxide 
take-up rate, mcc^ 1 are related to the biomass production 
rate, m bio ,by 

rho 2 = 1.0667 m bio 
m C o 2 =*-4m bi0 


The power requirement is determined using the lighting 
level recorded in table 2 measured as Photosynthetic 
Photon Flux (PPF), the radiation given off by the lights 
in the wavelength band useful for photosynthesis; the 
photoperiod, which is the hours each day that light is 
supplied to the plants; and the lighting system efficiency, 
Tj, which is the end-to-end efficiency of the lighting 
system. The equation used to calculate the required 
power, Preq> is 

PPF * photoperiod 

Preq = « 


We have not incorporated optimum lighting levels for 
crops in this study. For a legitimate application of these 
techniques, data representing crops at optimal conditions 
must be entered into table 2 or a more sophisticated, 
compensating model must be devised. Alternatively, 
multiple entries of the same crop species could be made 
with parameters reflecting the crop’s performance when 
optimized for biomass production, transpiration, power, or 
volume conservation. 

The final calculation is the total cost. This is the sum of 
the power and volume (height times lm^) times their 
corresponding penalties as given in table 1. 
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Total cost = (Preq) (Power penalty) 

+ (Volume) (Volume penalty) 

Note that this function could be made more sophisticated 
by weighting the relative importance of the two require- 
ments or by adding mass as an additional cost. Whereas 
mass penalties could easily be incorporated into the cost 
function, current understanding of mechanically opti- 
mized plant chamber mass is sufficiently limited to defer 
incorporating mass penalties into this study. 


A feature of the decision-tree tool is the accessibility of 
cost values for all crop mixes. This allows the designer to 
investigate the cost of nonoptimal solutions, which might 
have more appeal than the optimal solution for qualitative 
reasons. For the example shown in figure 2, increasing the 
crop variety by selecting the most optimal four-crop 
solution over the three-crop solution increases the cost by 
18%. Also, if a designer preferred wheat over potatoes in 
the optimal crop mix of lettuce, potatoes, and Y (crops 1, 
2, and 4), the decision tree shows an increase in cost of 
14% for lettuce, wheat, and Y (crops 1 , 3, and 4). 


Decision Analysis Method 

Decision analysis methods described in reference 14 
provide a tool for making decisions based on a single 
principal value (in our case, we have chosen to express 
everything in terms of cost). Most often, the tool used in 
decision analysis is a decision tree, where all possible 
outcomes and all possible paths to these outcomes are 
diagramed. Many decision-tree analyses also have 
expected values or probabilities attached to each branch 
stemming from a decision or node. In our case, we are 
merely minimizing cost at each decision node, with equal 
probability that any particular pathway will be followed. 
The decision tree (fig. 2) is constructed such that the 
initial decision determines which crop selection is the 
cheapest for all permutations of solutions having the same 
number of crops. A second round of decision-making is 
then done to determine the most economical number of 
crops one could use. 

Computation of the cost values is shown in table 4. The 
assumption is made that an equal area is allotted to each 
crop in a crop mix. Thus the generation/dissipation rates 
of the plant products, power, volume, and cost are 
averages of the individual crops in a crop mix computed 
on a per-square-meter basis. The penalties incurred by this 
assumption are shown through the comparison of results 
with the results from the design optimization method, 
where this assumption is not required. The required 
generation/dissipation rates (calculated from the human 
requirements and the degree of support specified in 
table 1) are divided by the productivity of each crop mix 
to obtain the scaling factor (planting area in square 
meters) required to meet the specified level of life 
support. The largest scaling factor encountered for a 
particular crop mix is multiplied by the cost of the crop on 
a per-square-meter basis to obtain the cost penalty entered 
in the decision tree. 


Design Optimization Method 

An alternative to using the decision analysis method 
described above is to take a design optimization approach 
(refs. 15 and 16). One can then minimize the cost function 
while removing the assumption of equal crop growth 
areas. For example, we use a linear programming 
approach to solve the constrained optimization problem 


N 


min< 

area 


(cost)j (area)j > 
i=l 


Subject to the following constraints 
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protein lower bound 
N 

protein produced j (area) . 

i-l i 

< protein upper bound 


carbohydrate lower bound 
N 

V" 1 f carbohydrate produced ^ 

~ Zj v area j ^ area ^ 

i-l i 

< carbohydrate upper bound 


lipid lower bound 

4 c ip - d ^ ) (area ,, 

i=l i 

< lipid upper bound 


where i identifies the crop species and N is the total 
number of crops being considered. The objective is then 
to identify the crop mix that minimizes cost while meeting 
certain requirements for air and water regeneration as well 
as food production. One could easily modify the above 
formulation to include additional requirements, such as 
vitamin and mineral nutritional requirements, or to 
include additional constraints, such as physical constraints 
on the crop-growth areas due to rack-size limitations or 
edge-effect considerations. 

This constrained optimization problem can be solved 
using standard linear programming techniques. We used 
a SIMPLEX algorithm, coded in FORTRAN on a 
MicroVax 3200 computer. The algorithm first identifies 
whether a feasible solution exists, then solves for the 
optimum solution and determines whether or not the 
solution is degenerate (i.e., an infinite number of solutions 
exist). 

The optimization method has computational advantages, 
especially when a large number of crop species or 
nutritional requirements are being considered. The 
decision analysis method requires an exhaustive search 
since the cost of each possible solution must be calcu- 
lated, whereas the optimization method uses search 
directions to quickly find the optimum solution. Also, 
standard methods exist for performing parameter 


sensitivity analyses for the linear programming 
formulation. Such analyses would be very useful for 
performing “what if’ studies to investigate the effects of 
changing costs or productivities of the various crop 
species. 

Results 

Table 5 shows the decision analysis output for the 
baseline case outlined in tables 1-3. For this case we have 
specified 100% of the requirements for oxygen, carbon 
dioxide^ and carbohydrates to be fulfilled, as well as 40% 
of those for water, 25% of lipids, and 25% of protein. 
Results show three crops being selected (lettuce, potatoes, 
and Y) as the optimum mix, with a total area of 3.0 m 2 
and total cost of 2.33. Results also show that in order to 
supply 100% of the carbohydrate requirement, oxygen, 
carbon dioxide, water, and protein are oversupplied, and 
more carbon dioxide is taken up than is necessary. In 
design optimization terminology, carbohydrates are the 
active constraint. 

Table 6 shows the results of the design optimization 
technique for the same base set of constraints. Recall that 
this technique is not limited to equal areas for each 
species in the crop mix selected. This is reflected in the 
results, which in this case show three crops being selected 
(lettuce, Y, and soybeans), with a total area of 2.5 m 2 at 
a cost of 2.16. Examining the results, we see that the 
carbohydrate, water, and lipid requirements are active 
constraints. 

The optimization technique can also be used to examine 
the optimal solution in the case where protein, carbo- 
hydrate, and lipid production all become active con- 
straints. This is accomplished by setting the upper and 
lower bounds on these variables all equal to 100% of the 
requirements, thus forcing the optimization to choose a 
crop selection which produces a nutritionally balanced 
diet (no overproduction of protein, carbohydrates, or 
lipids). In the baseline case, no feasible solution exists. 
However, if crop X is substituted for crop Y a solution 
exists such that exactly 100% of the protein, carbohydrate, 
and lipid requirements are satisfied. In this example, 
lettuce, potatoes, and crop X are selected, and carbon 
dioxide, oxygen, and water requirements are satisfied but 
are not active constraints. 

A study was performed by McDonnell Douglas using 
decision analysis to determine the optimum crop mix in a 
CELSS to fulfill nutritional requirements only (ref. 4). 
Their optimization is based on minimizing crop area (no 
penalties for power and volume consumption). In their 
study, it was found that the lowest total area results when 
no lipid production requirement is placed on the plants. In 
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the baseline case of our study, the lipid requirement to be 
satisfied by the plants is set to a low level (25%). For this 
case the lipid requirement is not an active constraint. 
However, when the requirements for the three nutritional 
categories are increased to a level of 100%, lipids become 
the active constraint, which is consistent with the 
McDonnell Douglas study. 

Sensitivities 

There are several ways in which sensitivities can be 
examined. Two approaches are demonstrated below. 

From the results of the baseline case, it was determined 
that carbohydrates were an active constraint. With this in 
mind, we might like to see how sensitive our answers are 
to the specified carbohydrate requirement. For example, 
what happens when the carbohydrate fulfillment level is 
lowered to 50%? The decision analysis results show a 
40% savings in cost and a 15% decrease in crop area 
using only soybeans and Y. With this new requirements 
specification, the active constraint shifted to water 
fulfillment. The design optimization method showed a 
38% savings in cost and a 21% decrease in crop area. The 
active constraints remain carbohydrate and water fulfill- 
ment. If these cases were based on real plants, hardware, 
and mission parameters, this sensitivity analysis would 
suggest that resupplying a portion of the carbohydrates 
from Earth might be preferred to requiring the plants to 
produce the entire requirement, depending on resupply 
costs. 

Another approach to analyzing sensitivities is to examine 
the effect on the results when plant performance data for a 
particular species are altered. As an example, the biomass 
production rate and transpiration rate were varied 
(doubled and halved) about nominal values. Table 7 lists 
the results of the two methods in each of these cases. This 
table shows that certain increases in plant performance 
parameters result in larger cost savings than others. This 
method could provide a tool for steering plant research, as 
well as addressing the balance between transpiration rate 
and plant biomass production. 

Conclusions 

1. The design optimization method inherently selects 
the most optimal crop mix of the two methods considered 
since it does not require all of the species selected to have 
equal crop areas. The decision analysis method results are 
usually within 10% of the cost and total crop area of the 
design optimization results. 

2. The advantages and disadvantages of both methods 
imply that the two methods are ideal companions to each 


other. The decision analysis method allows an interactive 
atmosphere with the user, facilitating experimentation 
with design specifications. A decision tree displays the 
cost of all possible crop mixes for the CELSS designers, 
allowing them to evaluate at a glance the cost of 
additional variety or preferred (more appetizing) crop 
species. 

The design optimization method has computational 
advantages, especially when the number of crops being 
considered or the number of nutritional (or other) 
constraints becomes large. It also ensures that the last 
10% of cost and area savings will be provided by the 
specified crop mix, and it allows the user to further 
constrain the optimization problem as needed. A more 
rigorous sensitivity analysis approach could also be 
developed using this method. 

3. There are limitations to this analysis. A significant 
limitation of the decision analysis method is the assump- 
tion of an equal area for each species in a crop mix. Both 
methods are limited because they do not allow the plant 
performance parameters to be simultaneously optimized 
(the transpiration rate/biomass generation rate trade-off, 
as well as other variables). No attempt is made to account 
for the compatibility of crop species if a common air 
space or nutrient delivery system is planned. Also, the 
power and volume of the processors required to support a 
plant system, including environment control, are not 
accounted for. 
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Table 1. CELSS specifications input 


Species to choose between (up to 5 ) 


LETTUCE 


POTATOES 


WHEAT 


Y 


SOYBEANS 


Minumum number of species 
required 1 


Portion of life support requirements 

fullfilled bv plants 

% Oxygen 

100 

% Carbon Dioxide 

100 

% Water 

40 

% Lipids 

25 

% Carbohydrates 

100 

% Protein 

25 


Penalties 

Cost/kW power 
Cost/m^ volume 


System parameters 


Lighting system 
efficiency 


1 Life support reauirements for 1 person 1 

Oxygen [Kg/day] 

0.84 

Carbon Dioxide [Kg/day] 

LOO 

Water [Kg/day] 

31.26 

Lipids [gm/day] 

105 

Carbohydrates [gm/day] 

633 

Protein [gm/day] 

76 































Table 2. Plant performance data for model input 


Height 

required 

Imi 

0.3 

0.11 

0.2 

0.6 

0.2 

0.3 

Lighting 

photoperiod 

[hr/dayl 

18 

20 

12 

20 

20 

20 

PPF 

level 

[umoi/m 2 s] 

485 

725 

532 

534 

400 

600 

Protein 
[% Tot. Bio.] 

22 

10 

40 

14 

5 

2 

Carbohydrates 
[% Tot. Bio.] 

58 

85 

20 

82 

65 

87 

Lipids 

[% Tot. Bio.] 

5 

0.5 

35 

2 

15 

3 

Biomass 

edible 

Fraction 

0.8 

0.8 

0.4 

0.5 

0.7 

0.5 

Biomass 
growth rate 
[Kg/m 2 day] 

0.6 

0.136 

0.3 

0.4 

0.3 

0.7 

Trans, rate 
[Kg/m 2 day] 

1.02 

6.41 

7.7 

9.98 

3 

5 

Plant species 

LETTUCE 

POTATOES 

SOYBEANS 

WHEAT 

X 

Y 


C3 

-a 

GJ 

O 

c 

C3 

E 


0> 

Oh 


PO 

S 

C3 

H 


Total 

cost* 

[$/m 2 ] 

0.7666 

0.6675 

1.3013 

0.8784 

0.5303 

Power 

tequired 

IkW/m 2 ] 

0.207 

0.344 

0.254 

0.285 

0.152 

o 2 

produced* 

[Kg/m 2 ] 

0.64 

0.15 

0.43 

0.64 

0.32 

* , 

Oh(n 

o 3 e 

u 1 3* 

03 > — i 

0.88 

0.20 

0.59 

0.88 

0.44 

Protein 

produced* 

[gm/m 2 ] 

105.6 

10.88 

28 

6 

48 

Carbohydrates 

produced* 

[gm/m 2 ] 

278.4 

92.48 

164 

261 

24 

Lipids 

produced* 

[gm/m2] 

24 

0.544 

4 

9 

42 

C/3 (— i 

V3 * CN 

i -s b 

1 1 3 1 

0.6 

0.136 

0.4 

0.6 

0.3 

•E rs -1 

D. * C 
c/3 0) c 

E 03 

2 ^ s/ 
H * 

1.02 

6.41 

9.98 

5 

7.7 

Plant species 

LETTUCE 

POTATOES 

WHEAT 

Y 

SOYBEANS 


03 
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Table 4 . Crop mix cost calculation 
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Table 5. Decision analysis baseline results 


Total crop area [m 2 ] 

3.0 

Number of plant species 

3 

Area planted per specie [m 2 ] 

1.0 

Total power required [kW] 

0.84 

Total volume required [m 2 ] 

1.40 

Total cost of plants 

2,33 


Species selected 


1 

LETTUCE 


2 

POTATOES 


3 

Y 


4 



5 




Life support requirements 
bv outout crop mix 

% Oxygen 

171 

% Carbon dioxide 

231 

% Water 

74 

% Lipids 

58 

% Carbohydrates 

100 

% Protein 

202 


Table 6. Design optimization baseline results 


Total cost 

= 

2.1595 



Total area 

= 

2.5004 

m 2 


Area of crop #1 (Lettuce) 

= 

0.0584 

m2 


Area of crop #2 (Potatoes) 

= 

0.0000 

m 2 


Area of crop #3 (Wheat) 

= 

0.0000 

m 2 


Area of crop #4 (Y) 

= 

2.3549 

m 2 


Area of crop #5 (Soybeans) 

= 

0.0870 

m 2 


C0 2 removed 

— 

2.1621 

kg/day 

(216.21 % of human reqt.) 

0 2 produced 

- 

1.5724 

kg/day 

(187.19 % of human reqt.) 

H 2 0 processed 

= 

12.5040 

kg/day 

( 40.00 % of human reqt.) 

Protein produced 

- 

24.48 

g/day 

( 32.21 % of human reqt.) 

Carbohydrate produced 

= 

633.00 

g/day 

(100.00 % of human reqt.) 

Lipid produced 

— 

26.25 

g/day 

( 25.00 % of human reqt.) 




Transpiration rate 

[kg/m 2 dl 


Crop 

Y 



2.5 

5.0 

(nominal) 

10.0 


0.35 

DA 

cost: 2.45 
area: 3.4 m 2 
species: L,P 

DA 

cost: 2.45 
area: 3.4 m 2 
species: L,P 

DA 

cost: 2.42 
area: 2.9 m 2 
species: L,Y 



OPT 

cost: 2.36 
area: 3.50 m 2 
species: L,S 

OPT 

cost: 2.36 
area: 3.50 m 2 
species: L,S 

OPT 

cost: 2.24 
area: 2.76 m 2 
species: L,Y 

Biomass 

generation 

rate 

[kg/m 2 d] 

0.70 

(nominal) 

DA 

cost: 2.42 
area: 3.3 m 2 
species: L,Y,S 

OPT 

cost: 2.28 
area: 2.98 m 2 
species: Y,S 

DA 

cost: 2.27 
area: 3.1 m 2 
species: L,Y,S 

OPT 

cost: 1.96 
area: 2.34 m 2 
species: Y,S 

DA 

cost: 1.87 
area: 2.3 m 2 
species: L,Y 

OPT 

cost: 1.79 
area: 2. 17 m 2 
species: L,Y 


1.40 

DA 

cost: 1.73 
area: 2.5 m 2 
species: Y,S 

DA 

cost: 1.41 
area: 2.0 m 2 
species: Y,S 

DA 

cost: 1.19 
area: 1.4 m 2 
species: Y 



OPT 

cost: 1.56 
area: 2.29 m 2 
species: Y,S 

OPT 

cost: 1.40 
area: 1.97 m 2 
species: Y,S 

OPT 

cost: 1.06 
area: 1.32 m 2 
species: Y,S 


LEGEND : 


P = 

potato 

DA = results from decision analysis 

L = 

lettuce 

OPT = results from design optimization 

W = 

wheat 


S = 

soybeans 


Y = 

Y 




















Figure 1. Conceptual design for a Controlled Ecological Life Support System. 
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Figure 2. Decision Tree. 
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